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ABSTRACT 

High-velocity features in Type la supernova spectra provide a way to probe 
the outer layers of these explosions. The maximum- light spectra of the unique 
Type la supernova 2000cx exhibit interesting Ca II features with high-velocity 
components. The infrared triplet absorption is quadruply notched, while the 
H&K absorption is wide and flat. Stimulated by a three-dimensional interpreta- 
tion of similar Ca II features in another Type la supernova (SN 2001el, Kasen 
et al. 2003), we present alternative spherically symmetric and three-dimensional 
ejecta models to fit the high-velocity {v > 16,000 km s^^) Ca II features of SN 
2000cx. We also present simple estimates of the high-velocity ejecta mass for a 
few trial compositions and discuss their implications for explosion modelling. 

Subject headings: supernovae: individual (SN 2000cx) — radiative transfer 



1. Introduction 

One-dimensional (ID) white dwarf explosion models like W7 (Nomoto et al. 1984) have 
been used to synthesize spectra in good agreement with those observed of Type la supernovae 
(SNe la; see Filippenko 1997, for a general review). But recent polarization observations of 
the SNe la 1997dt (Leonard et al. 2000b), 1999by (Howell et al. 2001), 2000cx (Leonard et 
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al. 2000a), and 2001el (Wang et al. 2003; Kasen et al. 2003) indicate that at least some SN 
la envelopes deviate significantly from spherical symmetry. 

The strongest case for such deviation appears in the Ca II lines of SN 2001el. This 
event exhibits an unusual Ca II infrared (IR) triplet in its flux spectrum about a week 
before maximum brightness. Two triplets are evident: one corresponds to photospheric- 
velocity (PV) material and the other to higher-velocity (HV) material. In the polarization 
spectrum, the HV feature coincides with significant intrinsic net polarization. Kasen et al. 
(2003) investigate a number of envelope models to account for the HV feature. Generally, 
they conclude that (1) incomplete covering of the photosphere and (2) some deviation from 
spherical symmetry are required to produce the HV features. Unfortunately, conditions did 
not permit simultaneous observation of the blue spectrum so a corresponding phenomenon 
in the Ca II H&K feature remains unconfirmed. 

Another unusual Ca II feature is found in the spectrum of the peculiar SN 2000cx (Li 
et al. 2001). Near maximum fight, its ostensible Ca II IR triplet possesses an interesting 
quadruply notched feature perhaps due to what Li et al. call "some unique distribution of 
Ca in the ejecta of SN 2000cx," extending to high velocity. The wavelength coverage of the 
near-maximum light spectra is excellent. Though polarization data at the same wavelengths 
are not available, simultaneous fits of the H&K and IR triplet features may help constrain 
at least the HV photospheric covering fraction (if indeed the HV material is not distributed 
in spherical symmetry). 

It remains unclear why clumpy HV ejecta could occur in SNe la. Detailed synthetic 

spectropolarimetry has yet to test explosion models free from the constraint of spherical 
symmetry (Khokhlov 2000; Reinecke et al. 2002; Gamezo et al. 2003). These explosion 
calculations have neither proceeded to the free-expansion phase nor provided abundance 
distributions required for detailed spectrum synthesis. On the other hand, methods for 
calculating synthetic spectra from three-dimensional (3D) models of SNe are in their infancy 
(Thomas et al. 2002; Kasen et al. 2003). Hence parameterized, direct analysis of observed 
SN spectra remains a powerful way of guiding explosion modellers to replicate geometrical 
phenomena in their models. 

In this article, the goal is to analyze the Ca II features in the spectrum of SN 2000cx at 
a single epoch to constrain the HV ejecta geometry. A separate multi-epoch, but exclusively 
ID analysis of this object is forthcoming (D. Branch et al., in preparation). That work 
addresses such orthogonal problems as the exhaustive identification of PV features and the 
unusual color evolution of this SN. 

The outfine of the remainder of this article is as follows. In §2 we present the spectra 
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and identify the features of interest. In §3 we use the ID direct analysis code Synow and the 
analogous 3D code Brute to test a few geometrical distributions of Ca II line optical depth. 
We discuss the candidate models and estimate the HV ejecta mass in §4. Conclusions appear 
in §5. 

2. Spectra 

In Figure 1 are three spectra of SN 2000cx collected near maximum hght, originally 
presented by Li et al. (2001). The usual SN la absorptions from Si II (near 6150 A) and 
S II (near 5400 A) are accompanied by absorptions due to Fe III (near 4300 A and 5000 A). 
There is very little, if any, signature of Fe II. 

Between 7900 A and 8400 A is a series of four notches. Usually at this phase, only 
absorption from the two stronger lines of the triplet (AA8542, 8662 — often blended) are 
visible in this region, and these are the best candidates for the two redder notches. The 
two bluer ones are perfectly consistent with a Ca II IR triplet shifted to about 22,000 km 
s^^ toward the observer. The bluest and weakest line of the Ca II IR triplet (A8498) is 
approximately ten times weaker than A8542, and is likely heavily blended with that line. 
Other ions are unlikely to produce the notches; candidates such as O I are unconvincing due 
to the absence or weakness of concomitant lines in the spectra. 

The corresponding Ca II H&K absorption feature (3500 A to 3800 A) is wide and flat. A 
collection of currently unidentified narrow absorption features obliterates its emission peak. 

Figure 2 displays the Ca II features of SNe 2000cx and 1994D near maximum light in 
terms of velocity relative to the observer. The top half of each panel is the H&K feature 
relative to the gfj-weighted doublet wavelength (3945 A), while the bottom half is the IR 
triplet relative to its g-/- weighted wavelength (8579 A). 

The fact that these two features are blends complicates the issue of choosing a reference 
wavelength for producing Doppler-space plots. In a blend of P-Cygni lines, redder compo- 
nents dominate the shape of the aggregate profile if they are strong; they screen out the bluer 
fines. Nevertheless, the bluer fines must have some effect since they feed radiation to the 
redder lines. Given these ambiguities, we compromise by using the g^J-weighted wavelengths 
for reference. 

Figures 2a, 2b, and 2c are from SN 2000cx at days 2, 6, and 7 past maximum light, 
respectively. For comparison, the normal SN la 1994D Ca II features at day 3 past maximum 
appear in Figure 2d. Hatano et al. (1999a) use HV Ca II and Fe II to improve their synthetic 
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fits of SN 1994D. In that SN, the effect of the HV Ca II is much less pronounced than in SN 
2000cx. In fact, the IR triplet of SN 2000cx has a HV absorption with a depth relative to 
the continuum between that of SNe 1994D and 2001el. 

Note that the four notches in the SN 2000cx IR triplet absorption do not evolve ap- 
preciably in velocity space over time, and that the overall velocity ranges of both features 

match. There appears to be a case for a onc-to-onc correspondence between notches in the 
IR triplet and depressions in the H&K absorption, but the blended nature of both features 
makes absolute confirmation difficult. Furthermore, the mild fiuctuations visible in the H&K 
absorption could easily be due to some weak, narrow lines superposed on the Ca II feature. 

We henceforth designate the parts of the Ca II features forming above 16,000 km s~^ 
relative to the (7 /-weighted wavelengths of H&K and the IR triplet as the "HV features." 
Our fitting strategy will be to concentrate on replicating the velocity range of these features 
and some of their structure. Since the three SN 2000cx spectra in Figures 1 and 2 are all 
quite similar, we restrict the focus of the remainder of this paper to the spectrum obtained 
two days past maximum light. 

3. Fits 

The simplest conceptual model of a SN la between a few days to months after explosion 
consists of an homologously expanding envelope surrounding a continuum-emitting, optically 
thick core. The decay of freshly synthesized ^^Ni to ^^Fe at the center of the SN releases 
energy in the form of 7-rays. This energy diffuses outward and degrades to ultraviolet, 
visible, and infrared wavelengths via interactions with the matter in the SN atmosphere. 

This elementary model has been implemented in a variety of codes intended for the 
empirical analysis of SN spectra (Jeffery & Branch 1990; Fisher 2000; Mazzah & Lucy 1993). 
These codes all implement a sharply defined photosphere at some radius as a lower boundary 
for a line-forming envelope. They utihze the Sobolev approximation (Sobolev 1947; Castor 
1970; Rybicki & Hummer 1978) to construct radiation field estimates for given optical depth 
distributions. Optical depths may be parameterized (in velocity space and wavelength) or 
determined from a selected abundance and density model with or without a self-consistent 
temperature structure. 

These codes have proved useful for fitting observed spectra to constrain the structure 
of SN envelopes. This empirical process is called "direct" analysis to distinguish it from 
"detailed" analysis where the full radiative transfer and non-local thermodynamic equilib- 
rium (NLTE) rate equations are solved (Hauschildt & Baron 1999). Specifically, the purpose 
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of direct analysis is (1) to determine what species are present in a SN line-forming region 
and (2) to constrain the velocity space distribution of those species. The assumptions and 
approximations used generally restrict direct analysis to line phenomena, in particular to 
the Doppler shifts and overall shapes of absorption features. The results of direct analysis 
provide useful guidance to detailed spectrum modellers and to explosion modellers as well. 

Synow is a direct analysis code that has been used to fit many spectra of various types 
of SNe, e.g., SN la 1994D (Hatano ct al. 1999a), several SNe lb (Branch et al. 2002), SN Ic 
19941 (Millard et al. 1999), and SN 11 1999em (Baron et al. 2000). Synow rehes on spherical 
symmetry and parameterized line optical depth. Each selected ion is assigned an optical 
depth in a "reference" line, usually a strong optical line. All other line optical depths for 
the same ion are determined assuming Boltzmann excitation of the levels at some specified 
temperature. The optical depths are scaled as a function of radius to fall off exponentially 
or according to a power law. The full details of Synow are described by Fisher (2000). 

We have recently developed a code similar to Synow in its degree of parameterization, but 
free from the constraint of spherical symmetry. This code. Brute, is based on Monte Carlo 
techniques presented in a series of papers by Lucy and Mazzah (e.g., Mazzali & Lucy 1993), 
and will be discussed in detail elsewhere (R. C. Thomas, in preparation). In Brute, spatial 

parameterization of line optical depth is managed through a 3D template. To establish 
the radiation field estimates in all transitions throughout the envelope, Monte Carlo energy 
packets are emitted from the core and their scattering histories arc followed. The required 
memory for Brute is much larger than that for Synow, and Brute also lacks the speed of its 
ID counterpart. 

3.1. Synow Fits 

Here we investigate a few spherically symmetric distributions of Ca II optical depth. 
The assumption of homology in the SN atmosphere {v oc r) permits us to parameterize 
Sobolcv optical depth in terms of velocity v relative to the explosion center. In spherical 
symmetry, we may designate domains Vmin < v < Vmax within which we define reference line 
Sobolev optical depth according to the rule 

Tref{v) = TrefiVmin) Q'^V[{Vmin " v) / V^] {Vmin <V < Vmax): (l) 

where Vg is an e-folding length. If Vmin is greater than the velocity at the photosphere Vph, 
we say the optical depth profile is detached. Outside of the domain we may set the reference 
optical depth to zero, or use another rule for the same reference line to set up a superposition 
of optical depth profiles. Other line optical depths for the same ion are assigned assuming 
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excitation temperature Texc- The velocity at the photosphere in all the presented Synow fits 
is Vph — 12, 500 km s~^. 

Synow uses a blackbody-emitting photosphere, clearly insufficient to account for all real 
continuum processes at work in the formation of a SN spectrum. This limits the range over 
which the synthetic continuum level can be made consistent with that observed, making fits 
of entire spectra extending from 3000 A to 10,000 A problematic. We adopt a piecewise 
approach and choose a convenient blackbody temperature applicable for the blue part of the 
spectrum. Fitting only the major features blueward of 6000 A (excluding lines of Ca II), we 
find that a blackbody temperature T55 = 12, 000 K yields a decent fit to this region. Since 
increasing T55 from this value changes the continuum slope in the neighborhood of the IR 
triplet only slightly, we choose to merely scale the synthetic spectra down in order to fit that 
feature. 

Except for the optical depth parameters of Ca II, the other parameters {vph, Tbb, and 
optical depths of ions fisted in Table 1) stay fixed from fit to fit. Since the Ca II IR triplet 
may blend with some lines of O I, a small amount of optical depth for that ion is included in 
the fit, but its presence has little impact on the results. The parameters used for the various 
Ca II optical depths are listed in Table 2. 

3.1.1. One- Component Fits 

Here we present three illustrative fits of the Ca II features using only one spherically 
symmetric velocity component in each. Figure 3 shows an example of a fit to the entire 
spectrum range. The observed features between 4500 A and 6000 A are fit rather well 
assuming the parameters for Fe III, Si II, and S II given in Table 1. 

In Figure 4 are close-ups of the H&K feature and IR triplet fits using model IDl from 
Table 2. This single-shell model with nearly constant optical depth as a function of radius 
can reproduce the broad velocity extent of the H&K feature. Yet it cannot reproduce any of 
the structure present in the IR triplet, so next we consider independent fits to the PV and 
HV features to investigate this structure. 

The appearance of the two strongest IR triplet fines (Ca II AA8542, 8662) in the PV 
feature as two distinct notches presents an interesting problem under spherical symmetry. 
In general, without assuming Vmin > Vph or imposing a finite Vmax, generating such a feature 
in spherical symmetry is only possible if Vph ^ Vgep, where Vgep is the velocity separation of 
the two lines. To fit the Si II and Fe III features, Vph = 12, 500 km s^^ is used, but Vgep for 
the two strongest IR triplet lines is only about 4000 km s~^ and both notches have minimum 
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wavelengths consistent with a 12,500 km s ^ blueshift. 

An alternative is to impose a finite v^ax or small that prevents blending of the two 
absorptions into one. But this has the undesirable effect of making absorptions with fiat 
bottoms which (when combined together) generate a Ca II A8542 feature shallower than the 
redder line, even though the former has a higher oscillator strength. For the moment, we 
allow the features to blend into one absorption, using the parameters listed for fit IDIPV in 
Table 2 which provide a satisfactory fit to the PV features in Figure 5. 

Fitting the HV notches with a shell of Ca II optical depth (fit IDIHV in Figure 6) 
is considerably less problematic than in the case of the PV feature. Using a small Ve or 
imposing a finite v„iax prevents blending that would otherwise unite the two bluer notches 
of the IR triplet feature. 

3.1.2. Two- Component Fit 

Combining the one-component PV and HV fits described above yields fit 1D2 (Figure 
7). This fit is satisfactory for the IR triplet, but its major deficiency is that the peak between 
the two blue notches is higher than observed. Adjusting v„iax from 25,000 km s~^ to higher 
velocity permits the redder feature to weaken this peak, but extends the bluest synthetic 
notch further to the blue than is desired. On the other hand, the synthetic H&K feature 
appears to need some higher-velocity material to extend its blue edge. 

3.1.3. Three- Component Fit 

Using three velocity components of Ca II optical depth allows us to fit every notch in 
the IR triplet absorption by simply overlapping pairs of notches formed by each component. 
The result of this model (1D3) is displayed in Figure 8. For the IR triplet, this appears to 
be the "best fit" among the ID models, but the H&K feature is too strong in its bluest part. 
Imposing a finite Vmax to counteract this HV tail only deepens the synthetic H&K absorption 
since it removes material that scatters light from the emission lobes of the envelope. 

3.2. Brute Fit 

The assumptions used in Brute are roughly the same as those in Synow except for the 
freedom to choose spatial optical depth distributions without spherical symmetry. Again the 
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same blackbody emitting, sharp photosphere is assumed, resulting in the same continuum- 
level problem described before, and we use the same remedy here. The technique used for 
radiative transfer in Brute is that of Monte Carlo, a difference from Synow, but both codes 
take into account multiple scattering in the calculation of the source functions of the lines. 

The spatial parameterization used for the Ca II optical depth consists of two compo- 
nents. One is spherical (PV material) and the other is not (HV material). The idea is to 
engineer a simple 3D distribution of optical depth that yields a synthetic spectrum consistent 
with observation, and fits at least as well as the ID model. 

For the HV material in the 3D fit, we adopt a simple geometry consisting of a circular 
cylinder of radius 8000 km s~^ coaxial with the line of sight to the center of the photosphere. 
The reference optical depth for Ca II inside the cylinder is assigned according to the rule 



According to this rule, the optical depth in the cylinder reaches a maximum value of T^ef — 63 
at the plane — 22,400 km s"^. Outside the cyhnder, the optical depth for the Ca II 
reference line is prescribed by eq. (1) with Tref{vmin) — 7, Ve — 3000 km s~^, and Vmin — 
13, 000 km s~^. A plot of reference optical depth along the line of sight is presented in Figure 
9. The reference optical depths of the other ions are parameterized exactly as in the Synow 
fits and produce the same features. 

The Ca II features resulting from this 3D parameterization are shown in Figure 10. The 
fit to the HV IR triplet is rather good. The synthetic HV H&K is more consistent with the 
observation than in the 1D2 or 1D3 parameterizations. Since the PV Ca II optical depth is 
distributed much the same as in the 1D2 case, the weak observed peak between the redder 
notches cannot be reproduced. In the ID case, this could only be accomplished by adding 
a third shell of optical depth and tuning its position. We note that in the 3D case, if the 
PV material is deployed in a nonspherical manner (such that parts of the photosphere shine 
through), then the peak between the two PV notches can be reproduced quite easily. 




0.0225v^ - 441.0 
-0.0315v^ + 768.6 



19, 600 <v^< 22, 400 km s 
22,400 < f;;, < 24,400 km s 



-1 



-1 



(2) 



4. 



Discussion 



4.1. HV Ejecta Geometry from the Fits 



Reconstructing SN envelope structure from spectra is an ill-posed inverse problem. The 
problem becomes even more difficult when the assumption of spherical symmetry is relaxed; 
the presented 3D solution represents but one of many possible solutions. 
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In fact, to first order, any HV distribution of Ca II optical depth which yields the 
same photospheric covering factor as a function of Vz will produce the same absorption 
features. The most promising avenue for breaking this degeneracy is through synthetic 
spcctropolarimetry of high-quality data. The Stokes polarization parameters Q and U for 
wavelengths covering the Ca II features are required. Such Q-U plots are the means by which 
the HV ejecta geometry of SN 2001el is constrained (Kasen et al. 2003). 

Both that study and the present one focus on the absorption features of Ca II, and ignore 
the emission features. Generally, emission phenomena are of limited utility in assessing HV 
ejecta geometry. In spherical symmetry, a HV shell gives rise to a weak, fiat-topped emission 
feature barely detectable in flux spectra. On the other hand, if the HV ejecta are organized 
into isolated regions of high optical depth, those not covering the photosphere along the 
line of sight will only contribute individual emission bumps to the spectrum. These weak 
features are even harder to detect if the PV ejecta are spherically symmetric and contribute 
significantly to the emission feature. 

On the practical side, the peculiarities of the spectrum of SN 2000cx make it impossible 
to use the Ca II emission features for geometrical analysis. The H&K emission feature is 
disrupted by a sequence of narrow lines, while the IR triplet emission is plagued by fringing 
in the CCD. 

These factors (the inherent degeneracies of the problem, the lack of corroborating po- 
larization data, the limited utility of emission features) prevent us from making all but the 
most conservative statements about the HV ejecta geometry of SN 2000cx. The similarity 
of the HV IR triplet of SN 2000cx to that of SN 2001el is suggestive, and its geometrical 
interpretation in the latter event motivates our consideration of a 3D model for SN 2000cx. 

Prom a purely empirical standpoint, none of the presented fits stands out as the "best" 
overall. Any cosmetic details missed by the synthetic spectra could be tuned away by making 
minor adjustments to the models. However, the main goal is to simultaneously fit the HV 
H&K and IR triplet, and the 3D fit seems to do a better job. 

Regardless of the details of the true geometry of the HV Ca II in SN 2000cx, the presence 
of this feature in only some SNe la evokes a series of questions to address in the future as 
high-quality data become available and 3D spectrum analysis techniques improve. 

(1) How frequent is HV ejecta clumping in SNe la? All SNe la might possess HV 
deviation from spherical symmetry which could escape detection if optically thick parts did 
not obscure the photosphere. SNe 2001cl and 2000cx might be examples where fortuitous 
orientation of HV ejecta permitted detection. This explanation is bolstered by the otherwise 
normal appearance of the SN 2001el spectrum. However, the peculiarities of SN 2000cx (the 
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Fe III lines in its spectrum, its unusual light curve and color evolution; Li et al. 2001) make 
this hypothesis somewhat problematic. 

(2) What imphcation does the presence of clumpy HV ejecta have for explosion models? 
If some SNe la have nonspherical HV ejecta while others do not, it might imply that more 
than one progenitor model is needed to explain the SN la phenomenon. More conservatively, 
the dumpiness of the HV ejecta might be a consequence of a so-called second parameter 
(Branch 2001). 

(3) Does the existence of nonspherical SNe la influence the prospects for precision cos- 
mology? The most obvious spectroscopic effects would be due to the apphcation of template- 
derived K-corrections (Nugent et al. 2002) from spherical SNe to nonspherical ones. The 
effect would be minor when the spectroscopic outcome of nonsphericity is isolated to a few 
lines (e.g. Si II). On the other hand, line-blanketing ions (e.g. Fe II) will greatly influence 
the photometry along some lines of sight more than others. 

Ca II lines are notoriously excitable over a large range of temperature and density. This 
makes it possible to place limits on the HV ejecta mass, and in turn to constrain future 
white dwarf explosion models. 



Here we make some simple estimates of the HV ejecta mass as constrained by the Ca II 
features visible in the SN 2000cx spectrum at two days after maximum light. We make 
two pairs of estimates, each pair consisting of a ID and 3D measurement. The first pair 
of estimates is made for the extreme case of a purely calcium HV ejecta composition. The 
second pair is based on a C/O-rich composition, similar to the outer, unburned layers of a 
white dwarf. The details of this composition are as listed in the SN ion signatures atlas of 
Hatano et al. (1999b). 

For each pair of estimates, we adopt the 1D2 and 3D parameterization domains for the 
volume used in calculating the HV ejecta mass: 



These presented mass estimates are not intended to be exact; in all cases the assumption 
of thermal equihbrium (TE) is employed. More rigorous constraints on the HV ejecta mass 
require detailed modelling of the physical conditions in the HV material, and these will be 
the focus of much future work. The estimates here and in Kasen et al. (2003) are a starting 
point for that effort. 



4.2. HV Ejecta Mass Estimates 




3.9 X 10^^ cm^, 
5.0 X 10^5 cm3. 



(3) 
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4-2.1. Pure Calcium Composition 

In the Sobolcv approximation, line optical depth for a transition / — > m in a SN envelope 
is given by (Jeffery & Branch 1990) 



7re^ 



fiuXmnM-^-^), (4) 



meC \ QuUi 

where is the oscillator strength, A/^ is the transition wavelength, ni and are the lower 
and upper level occupation number densities, respectively, and t is the time since explosion. 
The other symbols have their usual meanings. For convenience, we rewrite this as 

= 2.292 X 10-^(^/),„AAnt,^^Pt^|^, (5) 

where we have neglected the correction factor for stimulated emission, replaced the constants 
with their numerical values, and imposed TE. The partition function Q{T) for the ion in 
question is included. 

If the HV ejecta consist of singly ionized calcium in the ground state, then the optical 
depth in the Ca II A3934 transition (the reference line) becomes roughly proportional to the 
mass density of the HV ejecta. In reahty, the gas consists of other ions of calcium and other 
species, so deriving a mass density assuming a pure Ca II composition provides an extreme 
lower limit to the HV mass. Substituting A^ = 3934, {gf )3934 — 1.3614, Ei — eV, and 
td — 20 into eq. (5) gives 

r3934 ^ 2.AQncaii/QiT). (6) 

If we assume that the temperature in the HV ejecta is less than 8000 K, we can use the 
simple approximation T3934 ~ ncaii since Q{T) varies slowly from 2 to 3 up to 8000 K. 

Now an average mass density for the HV ejecta as derived from the Ca II optical depth 
can be written 

(Phv) - {Aca/NAvo)nca - 6.66 X W'^^thv) gcm-3. (7) 

If we average the reference line optical depths over the 1D2 and 3D parameterization domains, 
then (rj/y) = 25.5 and (rf/y) = 31.5. Substituting these values into eq. (7) and multiplying 
by the volume in each case, 

Mhf > 3.3xlO-8Mo, 

Mffv ^ 5.2 X IO-^Mq. ^ ^ 

These mass estimates are the "rock bottom" numbers required for the observational signature 
of the Ca II A8542 hne. The deflagration model W7 (Nomoto et al. 1984; Branch et al. 1985) 
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suggests that the amount of material and the densities in this region should be much higher, 
on the order of 6.0 x 10~^ Mq in spherical symmetry above 20,000 km s~^. Including other 
species into the mass calculations improves the lower limit, but this requires a choice of a 
particular composition. 



4-2.2. C/O-Rich Composition 

A more realistic composition model than that of pure calcium can be used to produce 

a more meaningful constraint on the HV ejccta mass. Here we adopt the C/O-rich compo- 
sition from the ion signature atlas of Hatano et al. (1999b) as the candidate model. This 
composition is representative of the unburned outer layers of an exploded white dwarf. 

The procedure for constraining the HV mass density is quite simple. For a given tem- 
perature T and an atomic line, we iteratively solve the equation of state for the mass density 

Ptot where r = 1 in the line. Repeating this procedure for a range of T yields a curve in the 
ptot — T plane. Generally, the space above the curve represents densities for which the line 
optical depth exceeds unity. In principle, if a line has a spectroscopic signature, then r > 1 
in the line and the corresponding ptot — T curve gives a lower limit to the mass density. Of 
course, this procedure depends on the assumption of thermal equilibrium, so for some ions 
the curves are not representative. 

Contours of r = 1 for Ca II A3934 and A8542 are shown along with curves for reference 
lines of other ions in Figure 11. For the C/O-rich composition used, the absolute minimum 
mass for which Ca II A8542 has a spectroscopic signature is at p^{^ = 2.2 x 10"-*^^ g cm~^. 

The value of p^in depends on the choice of T, as does the relative strength of Ca II 
A8542 to A3934. However, since the optical depth of Ca II A3934 is always much greater than 
that of the IR line, it is difficult to constrain the exact value of T in the HV ejecta using this 
technique. As r in a line approaches oo, the line profile saturates and the dependence of the 
absorption feature depth on r disappears. Rather than resort to predictions from detailed 
models for a temperature in the HV ejecta (as was done for SN 200 lei by Kasen et al. 2003), 
we adopt the absolute minimum value given above for a conservative estimate of the HV 
ejecta masses: 

Mhf > 4.3xl0^3^e, 

Mfj^ > 5.5 X IO-^Mq. ^ ' 

These mass estimates and densities arc roughly consistent with the W7 estimate of 6.0 x 
1O~^M0 in spherical symmetry above 20,000 km s~^. This might suggest that the Ca II 
signature could arise from primordial material as was suggested for SN 1994D (Hatano et 
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al. 1999a), rather than from freshly synthesized material. 

Other ion signature curves are plotted in Figure 11 along with the two Ca II curves. 
Portions of those curves lying above the Ca II A8542 curve indicate candidate ions to look for 
in the HV material. Of particular interest is the Fe II curve: Hatano et al. (1999a) detected 
the presence of Fe II in SN 1994D at HV. Another interesting ion is Ti II, usually detected 
as a trough of lines in SN 1991bg-like events (Filippenko 1992). The notches just to the blue 
of the 4300 A Fe III feature could be attributed to HV Ti II. 

One might consider using the non-detection of features of certain ions to place an upper 
limit on the mass at HV. This is difficult since many of the lines are blended with other, 
stronger features. For example, HV Sr II blends with the Ca II H&K feature. Lines from 
Fe II, especially if weak, blend with the Fe III features. Since we cannot be certain these 
lines are present, though weak and blended, we refrain from capping the HV mass estimate. 

4-2.3. Other Compositions 

Marietta et al. (2000) present interesting simulations of the effect that a white dwarf 
in a binary has on its companion when it explodes. Generally, they note that more evolved 
companions (with less tightly bound envelopes) are more vulnerable to losing a substantial 
fraction of their envelope during the explosion. The resulting ejecta distribution includes an 
evacuated cone in the ejecta behind the companion and some hydrogen with characteristic 
velocity on the order of 1000 km s~^. However, a small amount of stripped hydrogen (about 
IO^^Mq) could be carried up to velocities exceeding 15,000 km s~^. Could the HV Ca II 
lines be a signature of such material? 

Figure 12 is another ion signature plot, this time using the H-rich composition (Hatano 
et al. 1999b). The same Ca II signature curves are plotted, and the curves for both Ha and 
H/3 are included for comparison. The required minimum density for an optical depth of 1 
in Ca II A8542 is 45% higher than in the C/O-rich case, and the density for which the Ha 
optical depth becomes unity is only 4 times higher than that. This leaves a rather narrow 
density window in which to form the HV Ca II without also producing an Ha signature, if 
the Ha signature were to be produced under the assumption of TE. 

It is an observational fact, however, that Ha in SNe II is generally in net or total emission 
(see, for example, spectra of SN 1999em - Hamuy et al. 2001; Leonard et al. 2002; Elmhamdi 
et al. 2003), an effect which cannot be replicated by a direct analysis code which uses pure 
resonance scattering line source functions where line ffux is conserved except for occupation 
effects. However, pure resonance scattering is sufficient for simultaneous, satisfactory fits 
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of H/3 and H7 in these objects (e.g., Baron et al. 2000, where synthetic Ha is permitted to 
remain too deep compared to observations in Synow fits of SN 1999em). 

If the HV ejecta in SN 2000cx were a clump of HV H-rich material blocking the photo- 
sphere, then the usual Ho; emission phenomenon has an interesting consequence. Assuming 
pure resonance scattering for all the Balmer lines, this material would produce HV Ha and 

H/3 absorption features at the same velocities as the HV Ca II features. The blueshifts are 
such that the HV Ha absorption falls into the Si II feature near 6100 A. But the net emis- 
sion phenomenon will mitigate the strength of the Ha absorption and further conceal the 
hydrogen signature. We would expect the H/3 feature to remain in absorption, however. 

In Figure 13, we show a fit to the entire spectrum with and without hydrogen in the 
3D HV ejecta. The hydrogen reference line (Ha) optical depth used for this spectrum has a 
velocity-space profile similar to that described in eq. (2) except that the maximum optical 
depth at = 22,400 km s^^ is r^e/ = 4. As predicted, the assumption of pure resonance 
scattering produces a readily apparent modification to the Si II absorption feature which 
would be ameliorated by the net emission effect. Additionally, a previously unidentified 
feature at 4500 A is fit by the synthetic HV H/3. 

Note that this signature is only possible if the HV hydrogen is confined to a clump 
blocking the photosphere. If the HV hydrogen is placed into a HV shell, a strong emission 
feature centered at the Ha rest wavelength will appear. This is clearly not the case in the 
spectra of SN 2000cx. Other candidate ions for fitting this fine simply fail, as discussed in a 
ID direct analysis of this object (D. Branch et al., in preparation). This fit is included here 
since it is contingent upon the 3D distribution of hydrogen at HV that we describe. 

If the 4500 A feature is HV H/3, and we dispense with the assumption of TE for the 
hydrogen lines, then we see that the lower level of H/3 could potentially be overpopulated. 
This reduces the amount of HV ejecta required to produce an H/? feature by a factor equal to 
the departure coefficient of the lower level of this transition. This reduces the gap between 
the HV masses estimated in this work and those predicted by (Marietta et al. 2000). But 
can NLTE effects bridge this gap? Future detailed 3D NLTE calculations are required to 
answer this question. 

Returning to Figure 12, we note the other ions that could be detected if the HV ejecta 
were indeed H-rich. It is unfortunate that basically the same ions as in the C/O-rich case 
appear. The strongest discriminator between the two compositions is the Balmer hydrogen 
series. The majority of the hydrogen stripped from the companion in the Marietta et al. 
(2000) models is at very low velocity, so it seems that late-time spectra should reveal this 
material. But as Marietta et al. (2000) themselves point out, detection of the material is 
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problematic. Nebular Ha is severely blended with iron and cobalt lines. Resorting to the 
infrared introduces problems with atmospheric water lines (interfering with Pa), and other 
telluric and blending effects (interfering with P/3). 

However, He I does make an appearance (in TE) at high temperature. If this ion is 
nonthermally excited, it might produce an HV signature. While He I A5876 is not observed, 
nonthermal excitation could produce He I A10830. At HV this line could explain a feature 
previously attributed by Rudy et al. (2001) to Mg II. 

Another scenario for the production of a SN la is through recurrent novae (Starrfield 
2003) or supersoft x-ray sources (Hachisu et al. 1999). These progenitors involve a He- rich 
(and H-deficient) companion star. Minimum mass estimates for this case suggest HV ejecta 
masses about the same as in the H-rich case, assuming TE. But since helium stars are more 
tightly bound than those simulated by Marietta et al. (2000), mass stripping from a He- rich 
companion seems less likely an origin for the HV material. 

Arguments in favor of a circumstcllar origin fall short. The spectroscopic signature for 
that material is completely different from what is seen in the case of SN 2000cx. In that case, 
wc expect a narrow Ha emission spike, not unlike that found in SNe Iln (e.g., Filippenko 
1997). 

5. Conclusions 

We have presented several exploratory fits to the unique Ca II features of the unusual 
SN la 2000cx. A ID shell of material can account for the HV IR triplet feature, but has 
difficulty doing so for the corresponding H&K feature. In 3D, both HV features can be fit 
simultaneously with a chunk of material along the line of sight which partially covers the 
photosphere. 

Assuming a C/O-rich composition, mass estimates for the HV ejecta in both geome- 
tries discussed are consistent in the lower limit with the W7 model, suggesting a possible 
primordial origin. If the HV ejecta are H-rich, however, the mass required at HV to produce 
the observed Ca II features is somewhat higher. There is circumstantial evidence for the 
HV ejecta being H-rich, pointing to a signature of stripping from a companion star. This 
model is more favorable when NLTE effects are taken into account, but confirmation requires 
detailed 3D spectrum synthesis which is not currently feasible. 

Understanding the origin of the HV ejecta in SNe la (and also its frequency and physical 
conditions) is potentially quite important for future explosion models and for answering the 



-16- 



progenitor question. It is clear that better constraints on the HV Ca II phenomenon require 
high-quahty flux and polarization spectra at near-infrared and near-ultraviolet wavelengths. 
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Table 1. Synow Fit Parameters for Non-Calcium Ions 



'^refx^min) ^min "^max ^ exc 

{W km s-i) (lO^ km s'^) (lO^ km s"^) (lO^ K) 



Felll 1.5 12.5 oo 1.0 10.0 

Sill 3.5 12.5 oo 1.0 10.0 

SIX 1.5 12.5 oo 1.0 10.0 

1 0.2 12.5 oo 3.0 8.0 
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Table 2. Synow Ca II Fit Parameters 



Fit 


(Figure) 


Number of 






max 








Components 




(10^ km s"^) 


(10^ km s"^) 


(10^ km s"^) 


IDl 


(4) 


1 


1.4 


12.5 


30.0 


20.0 


IDIPV 


(5) 


1 


10.0 


12.5 


oo 


3.0 


IDIHV 


(6) 


1 


20.0 


2 1.0 




0.5 


1D2 


(7) 


2 


10.0 


12.5 


24.0 


3.0 








30.0 


24.0 


25.0 


3.0 


1D3 


(8) 


3 


16.0 


13.0 


15.1 


3.0 








7.0 


19.0 


23.5 


3.0 








12.0 


23.5 


oo 


3.0 
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Wavelength (A) 

Fig. 1. — Near-maximum light spectra of SN 2000cx (Li et al. 2001). The epoch relative to 
maximum light is listed above each spectrum. All spectra shown in this paper have been 
corrected for the redshift of the host galaxy, NGC 524, z = 0.0080. 
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Fig. 2. — Ca II features in SN 2000cx at days 2 (a), 6 (b), and 7 (c) after maximum light, 
and those in SN 1994D at day 3 (d) after maximum. The features are plotted in terms of 
velocity relative to the observer, using the /-weighted wavelengths of the Ca II H&K and 
IR triplet features. 
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Fig. 3. — Synow fit IDl (dotted line) of SN 2000cx (solid line) two days after maximum light 
with one Ca II component. Optical depths and excitation temperatures used for the other 
ions (Fe III, Si II, S II, and O I) are as hsted in Table 1. 
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Fig. 4. — Synow fit IDl (dotted line) of Ca II features in SN 2000cx (solid line) two days after 
maximum light. A single, nearly constant optical depth shell extending from the photosphere 
to 30,000 km s~^ is used. The velocity extent of Ca II H&K is approximately reproduced, 
but the synthetic IR triplet lacks any of the structure seen in the observed spectrum. 
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Fig. 5. — Synow fit IDIPV (dotted line) of Ca II features in SN 2000cx (solid line) two days 
after maximum light. A single, exponentially decreasing optical depth component just above 
the photosphere is used. 
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Fig. 6.— Synow fit IDIHV (dotted fine) of Ca II features in SN 2000cx (sofid fine) two days 
after maximum fight. A single, exponentiaUy decreasing optical depth component at 20,000 
km s~^ is used. An extremely small e-folding velocity of 500 km s~^ is needed to prevent 
blending. 
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Fig. 7. — Synow fit 1D2 (dotted line) of Ca II features in SN 2000cx (solid line) two days 
after maximum light. Two exponentially decreasing shells of Ca II optical depth are used. 
One extends from the photosphere to 24,000 km s~^ and the other from 24,000 km s~^ to 
25,000 km s-^ 
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Fig. 8. — Synow fit 1D3 (dotted line) of Ca II features in SN 2000cx (solid line) two days 
after maximum light. Three exponentially decreasing shells of Ca II optical depth are used 
to generate the peak between the two red notches of the triplet. The synthetic H&K feature 
is too strong compared to the observed one. 
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Fig. 9. — Ca II reference line optical depth along the —z axis in the 3D model. 
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Fig. 10. — Brute fit (dotted line) of Ca II features in SN 2000cx (solid line) two days after 
maximum light. The IR triplet is fit, and the blue edge of the synthetic H&K feature is 
weakened with respect to fit 1D3. 
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Fig. 11. — Contours of r = 1 for a variety of ions in the C/O-rich composition (Hatano et al. 
1999b), assuming TE. The Ca 11 signature curves are repeated in both panels for comparison 
with curves of other ions. 
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Fig. 12. — Contours of r = 1 for a variety of ions in the H-rich composition (Hatano et al. 
1999b), assuming TE. The Ca II, Ha, and H/3 signature curves are repeated in both panels 
for comparison with curves of other ions. 
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Fig. 13. — Brute fit to SN 2000cx spectrum at two days after maximum, with (solid line) 
and without (dotted line) H optical depth in the HV clump. Ha absorption is concealed by 
the Si II absorption at 6100 A. The efi^ect of Ho; could be completely erased by net emission 
eff'ects. 



